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PREFACE 

BAMIRAC,  the  Ballistic  Missile  Radiation  Analysts  Center,  Is  a  facility  of  the 
Institute  of  Science  and  Technology  at  The  University  of  Michigan.  It  Is  supported 
by  ARPA,  the  Advanced  Research  Projects  Agency,  under  Contract  SD-91,  as  a 
part  of  Project  DEFENDER  (research  on  and  defense  against  ballistic  missiles). 

BAMIRAC  functions  as  a  center  for  analysts  of  and  Information  on  ballistic 
missiles  and  space  vehicles,  with  the  objective  of  formulating  defense  measures 
against  them.  Analyses  of  experimental  and  theoretical  results,  conducted  in  three 
general  subdivisions  corresponding  to  the  launch,  midcourse,  and  reentry  phase 
of  ballistic  missile  flight,  are  evaluated  and  examined  for  correlations.  In  Its  opera¬ 
tion  as  a  technical  Information  center,  BAMIRAC  collects,  processes,  and  dissemi¬ 
nates  Information  concerning  electromagnetic  radiation  emanating  from  or  caused  by 
ICBM’s  or  IRBM's,from  reports  of  field  measurements  as  well  as  from  laboratory 
and  theoretical  results. 

In  addition  to  analyzing  the  results  from  external  sources,  BAMIRAC  conducts 
Its  own  theoretical  and  experimental  Investigations. 

The  results  of  the  work  within  BAMIRAC  are  combined  with  the  material  obtained 
from  outside  sources  to  develop  models  for  ballistic  missile  radiation.  The  Informa¬ 
tion  so  obtained  is  disseminated  in  technical  reports,  which  are  published  frequently. 
The  technical  content  of  the  library  holdings  is  made  available  to  organizations  whose 
representatives  have  proper  authorization ;  in  addition,  the  library  prepares  bibliog¬ 
raphies  on  request.  BAMIRAC  also  aids  ARPA  by  handling  various  technical  meet¬ 
ings  and  by  publishing  the  Proceedings  of  the  Anti-Missile  Research  Advisory  Council. 

BAMIRAC  Is  under  the  technical  direction  of  the  Infrared  Laboratory.  It  draws 
also,  however,  upon  the  capabilities  of  the  Computation  Department  of  the  Institute, 
and  upon  those  cf  the  Aircraft  Propulsion  Laboratory  of  the  Department  of  Aeronauti¬ 
cal  and  Astronautlcal  Engineering  and  of  other  departments  of  The  University  ot 
Michigan,  particularly  within  the  College  of  Engineering. 
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SYMBOLS,  CONSTANTS,  AND  CONVERSION  FACTORS 

"1 

A™  *  Einstein  transition  probability  from  state  mton,  sec'1 

a  =  Classical  semi-major  axis  of  electron  energy  state,  cm 

Sq.  a+  Defined  by  Equations  56  and  58 

C  Defined  by  Equation  42 

c  =  Velocity  of  light,  2,997929  x  10*®  cm/sec 

d  *=  Debye -Hilckel  (subscript) 

E  =  Term  value  of  state  m,  cm"* 

=  Energy  level  of  state  m,  ergs 

e  =  Electron  charge,  4.80266  x  10**®  esu 

g  *  Degeneracy  of  state  m,  2J  +  1 
m  m  _27 

h  =  Planck's  constant,  6.62517  x  10*  erg  sec 

-1  -3  -1 

I  =  Intensity  of  spectral  line,  erg  sec  cm  sterad 

J  =  Inner  quantum  number  of  state  m 

m  ,6  , 

k  =  Boltzmann's  constant,  1.38044  x  10  erg  deg 
-  1/1.4388  cm'*  deg"* 

Kp  «  Equilibrium  constant,  atm 

m  «  Energy  state 

•28 

m£  >  Electron  mass,  9.1083  x  10  g 

N  >  Original  mole  fraction  of  copper 

n  x  Total  particle  density 

-3 

n  x  Number  density  of  atoms  in  state  m,  cm 

m  .j 

ng  x  Electron  number  density,  cm 

n^  -  Number  density  of  atoms  in  Ionisation  state  q.  q  ■  0,  1,  2,  .  .  . 
p  ■  Pressure,  atm 

R  ■  Intensity  ratio 

T  ■  Temperature,  °K 

Vq  ■  Energy  to  Ionize  an  atom  from  q-th  to  (q  +  l)-at  ionization  state  In  the  plasma 
xB  ■  Mole  fraction  of  species  s 

Zq  x  Partition  function  for  atoms  in  the  q-th  ionisation  state 

a  *  (2irmek)^Vh^  (Section  6) 

"  -  Vzo 

v  ■  Frequency,  sec"1 

Ad  -  Debye  radius,  cm 

A  ■  Wavelength 
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ir  =  3.14159  .  .  . 

£  =  Effective  nuclear  charge  In  units  of  e 

Cu  I  =  Neutral  copper 
Cu  II  =  Singly  Ionized  copper 

Conversion  Factors  [1] 
lev  =  1.60206  x  10'12  erg 
=  8068.03  cm*1 
=  2.41814  x  1014  sec*1 


,  -1 

1  cm 

=  1.98618  x  10*16  erg 

=  12397.87  x  10*8  ev 

-1 

-27 

1  sec 

=  8.82517  x  10  erg 

=  4.13541  x  10*15  ev 

1  erg 

=  1.50940  x  l02fl  sec*1 

=  5.03479  x  1015  cm'1 

1  atm 

=  1.013248  x  106  dynes 

Constants  for  Copper  [7],  [8] 

For  Cu  I  X  5105.54 

g  Am  -  .051  x  108  sec"1 
m  n  * 

E  *  30783.7  cm 

m 

For  Cu  I  A  5153.24 

g  Am  *  4.7  x  10®  sec'1 
°m  n  , 

E  *  49935.2  cm* 

m 

Ionization  potential  for  Cu  I 

*  7.7  2  588  ev 

=  82317.2  cm'1 

Ionization  potential  Cu  n 

*  20.2907  ev 

■  183885.8  cm*1 


* 


PROPERTIES  OF  COPPER  PLASMA 
Calculation  of  tho  Partition  Function,  Spocios  Concentration, 
and  Spectral-Line  Intensities 

ABSTRACT 

The  self-consistent  partition  function;  ionization  potential  lowering;  electron, 
ion,  and  neutral- particle  species  concentrations;  and  the  absolute  spectral- line 
intensities  of  Cu  I  A  5106  and  Cu  I  A  5153  are  reviewed  and  tabulated  for  pressures 
between  0.03  and  10  atm,  for  original  mole  fractions  between  10*4  and  0.03  and 
for  temperatures  from  2000°  to  10,000°K,  for  copper  vapor  in  a  hot  inert  gas  at 
thermal  equilibrium.  Only  single  ionization  of  copper  is  considered  in  this  range 
of  conditions.  For  the  partition  function  cutoff,  all  electron  states  whose  classical 
semi-major  axis  exceeds  one-half  the  Debye  distance  are  suppressed.  The  results 
should  be  of  use  in  spectroscopic  plasma  diagnostics  and  other  situations  where 
the  trace  amounts  of  copper  usually  present  can  be  used  for  temperature  measure¬ 
ment. 


1 

INTRODUCTION 

For  the  determination  of  temperature  and  other  physical  properties  within  a  plasma,  spec¬ 
troscopic  methods  are  of  great  interest  and  value  [2,  3,  4] . 

In  this  report  the  formulas  necessary  to  calculate  species  concentration  and  spectral  line 
Intensities  are  reviewed  and  then  applied  to  copper  vapor  in  an  inert  gas.  These  calculations, 
while  straightforward  for  atomic  spectral  lines  in  optically  thin  plasmas,  are  rather  tedious, 
particularly  In  view  of  the  Iteration  necessary  even  for  this  rather  simplified  case;  therefore 
a  digital  computer  was  used. 

In  order  to  make  meaningful  estimates  of  the  required  sensitivity  of  photomultipliers, 
amplifiers,  and  other  light  detectors,  it  Is  necessary  to  have  quantitative  Information  on  the 
intensity  of  radiation  that  will  be  measured  for  the  range  of  temperatures,  pressures,  and  mole 
fractions  that  are  expected. 

The  intensities  of  the  spectral  lines  Cu  I  A  5106  and  Cu  I  A  5153  are  calculated  for  the 
idealized  case  of  copper  vapor  in  a  hot  inert  gas.  All  electrons  are  presumed  to  originate 
from  the  single  ionization  of  copper.  Double  ionization  of  copper  has  been  neglected  for  the 
range  of  conditions  covered  here.  This  assumption  gives  an  error  of  less  than  1/2%  at  10,000°K. 
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The  partition  functions  for  Cu  I  and  Cu  II,  ionization  potential  lowering,  Debye  radius,  and 
species  concentrations  are  also  tabulated.  This  allows  computation  of  other  copper  Bpectral  line 
Intensities  if  values  for  their  transition  probabilities  are  known. 

It  should  be  possible  to  use  the  calculations  to  help  determine  contamination  levels  of  copper 
electrode  material  Introduced  in  plasma  jets,  as  well  as  temperatures  within  the  plasma. 

The  internal  partition  function  tabulation  will  permit  calculation  of  the  contribution  of  the 
Internal  degrees  of  freedom  to  the  various  thermodynamic  functions. 


2 

FUNDAMENTAL  EQUATIONS 

2.1.  INTENSITY  OF  AN  ATOMIC  SPECTRAL  LINE 

The  radiant  power  from  a  spectral  line  emitted  in  all  directions  (4ir  sterad)  from  a  unit 
volume  of  gas,  arising  from  a  transition  between  two  energy  states  of  an  atom,  is  equal  to  the 
product  of  the  number  of  transitions  per  unit  time  between  the  two  Btates,  and  the  energy  change 
Involved  in  the  transition  between  the  two  states. 

The  number  of  transitions  per  second  between  an  upper  state  m  and  a  lower  state  n  is 

3 

given  by  the  product  of  the  number  of  atoms  per  cm  in  the  upper  state  m,  and  the  transition 
probability  per  second  from  state  m  to  state  n. 

-3 

If  the  number  density  of  atoms  In  the  upper  state  m  is  n^  cm  ,  and  the  transition  prob- 

ability  between  upper  state  m  and  lower  state  n  is  A  sec  ,  then  the  number  of  transitions 

n 

per  second  between  upper  state  m  and  lower  state  n  is 


.m  -3  -1 

n_A  cm  sec 
m  n 

The  energy  change  of  the  transition  between  an  upper  state  of  energy  E  and  a  lower  state 

m 

of  energy  Er  Is  given  by 

E  -  E  *  hv  ergs 
m  n 

where  v  is  the  frequency  of  emitted  radiation  and  h  is  Planck's  constant. 

3 

Therefore,  the  radtant  intensity  per  cm  of  gas  emitted  into  1  sterad  from  a  spectral  line 

is  [4] 

I  ■  -X  hvAmn  erg  sec"*  cm"^  sterad"*  (1) 

4  it  n  m 
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At  thermal  equilibrium,  the  number  density  of  particles  in  energy  state  1.  compared  to  the 
number  density  of  particles  in  energy  state  m,  Is  given  by  the  Boltzmann  energy  distribution 
as  [5,  6] 


n 


m 


-Et/kT 

e 

-E  /kT 
e  m 


where  g  Is  the  statistical  weight  or  degeneracy  of  energy  level  Ej,  and  k  Is  the  Boltzmann  con¬ 
stant.  Summing  over  all  energy  levelB  in  the  same  Ionization  state,  we  find  the  number  den 
slty  of  atoms  In  the  q-th  state  of  ionization  to  be  [4,  5] 


v_£  -vCT 

-E  /kT  V  g‘ 
m  i 


n  Z 
m  q 

■  E  /kT 


m 


-3 

cm 


(2) 


where 


Z 


q 


L 


-E./kT 

Kle 


(3) 


Is  the  Internal  partition  function  for  atoms  In  the  q-th  state  of  ionization.  Therefore, 


nn*m  6 

q  m 


-E  /kT 
m 


m 


cm 


-3 


(4) 


Substituting  Equation  (4)  Into  Equation  (1),  we  find  [4] 


m  'Em/W 
n  hv  g  A  e 
q  m  n 

Z  4rr 


-1  -3  .  .-1 

erg  sec  cm  sterad 


(5) 


Values  of  E  are  tabulated  In  Reference  7,  and  values  of  g_Am  can  be  found  In  Reference 
m  m  n 

8.  Other  compilations  and  bibliographies  of  transition  probabilities  can  be  found  In  References 

9,  10,  and  11. 
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2.2.  SAHA  EQUATION 

For  a  given  element,  the  number  (tensity  of  atoms  in  successive  ionization  states  at  ther¬ 
mal  equilibrium  is  given  by  the  Saha  Equation  [4,  12,  13): 


n  ,n  2Z  ,  (2irm  kT) 
q+1  e  _  q+1  e 

n_  =  Z_  ~5 


3/2 


-Vq/kT 


(6) 


where  n  =  number  density  of  atoms  in  the  q-th  ionization  state,  cm 

n 

ng  =  number  density  of  electrons,  cm‘J 

Zq  =  partition  function  for  atoms  in  q-th  ionization  state 

m  =  electron  mass 
e 

Vq  -  ionization  energy  necessary  to  ionize  an  atom  from  the  q-th  state  of  ionization  to 
the  (q  +  l)-st  ionization  state  within  the  plasma. 


Note  that  the  mass  of  the  atom  in  the  q-th  state  of  ionization,  niq,  has  been  taken  equal  to 

m  ,  and  Z  -  2. 
q+1  e 

Considering  single  ionization  of  neutral  atoms  (q  »  0),  Equation  6  becomes 


n  n  2Z  (2irm  kT)3/2  -Vn/kT 
+  e  +  e _  0 

~Z0  h3 


Within  the  plasma  we  have  charge  neutrality  [14] 


(7) 


n  » 
e 


L 


<in„ 


where  the  sum  is  taken  over  all  ionization  states,  q  =  0,  1,  2,  .  .  .  . 


(8) 


Equations  5,  6,  and  8  are  the  fundamental  equations  necessary  for  the  calculation  of  species 
concentration  and  spectral-line  intensities.  They  can  be  manipulated  in  several  ways  for  con¬ 
venience  in  a  particular  application. 


3 

COFFER  SPECTRAL  LINES,  Cu  I  X  5106  AND  Cw  I  k  S1S3 

The  Cu  I X  5106  and  Cu  I  X  5153  lines  are  of  interest  here  for  the  following  reasons: 

1.  Their  wavelengths  are  nearly  equal,  which  will  minimize  problems  of  spectral  sensi 
tlvlty  of  a  detector. 
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2.  Their  upper  energy  levels  have  relatively  large  differences,  so  that  the  population 
ratio  of  these  levels  Is  a  relatively  strong  function  of  temperature,  t  rough  the  Boltz¬ 
mann  factor. 

3.  Neither  line  Is  a  ground-state  transition;  thus  problems  of  self-reversal  of  the  central 
radiation  by  the  relatively  cold  outer  atoms  Is  avoided. 

4.  They  have  relatively  well-determined  transition  probabilities. 

5.  They  arise  from  neutral  copper  atoms  so  that  partition  functions  do  not  enter  Into  their 
relative  Intensity. 

6.  They  have  reasonable  Intensities  and  Intensity  ratios  for  temperatures  in  the  range  3000° 
to  10,000°K. 

The  Intensity  ratio  of  these  two  lines  Is  calculated  below.  Their  absolute  values  are  tab¬ 
ulated  In  the  appendix,  having  been  calculated  from  Equation  5. 


The  Intensity  ratio  of  two  lines,  1  and  2,  arising  from  an  atom  in  the  Bame  ionization  state 
(nQ1  -  nQ2,  ZQ1  -  ZQ2)  is,  from  Equation  5, 


,  m 


'l  .  'I  «ml  Anl  _ 

rmf - bT“1> 


(-Hr522) 


(») 


so  that 


'l  ._  V  («„An )l  E--E 


logy  log—*  log  - 

2  *  (vOs 


ml  m2 , 

IcT - l0ge 


Assuming  the  index  of  refraction  for  the  two  wavelengths  as  measured  in  air  is  nearly  con¬ 
stant,  then 

v\/v2  “  VA1 


If  E_  Is  measured  in  cm  ,  then 
m 


'l  ,_MT)l  >1  0M<8W%1  ' 


J  («mA»  k 


For  Cu  I  [7,  8]  (see  also  constants  in  Symbols  List), 


-  5153.24 
\2  -  5105.54 


(10) 
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4 

IONIZATION-POTENTIAL  LOWERING  AND  TRUNCATION 
OP  THE  PARTITION  FUNCTION  SUM 


Because  of  the  Interaction  of  an  atom  with  the  nearby  particles  in  a  plasma,  the  ionization 
potential  will  be  lowered  for  a  particle  inside  a  plasma.  Physically,  an  excited  atom  InBlde  a 
plasma  can  spontaneously  ionize  because  of  the  perturbation  of  its  neighbors  if  its  energy  state 
is  between  this  lowered  ionization  potential  and  the  ordinary  ionization  limit  of  an  unperturbed 
atom. 

Coupled  with  the  physical  process  described  above  is  the  observation  that  for  large  temper 
atures  the  exponential  factor  In  the  internal  partition  function  Bum  (Equation  3) 


Z 


gKe 


-EK/kT 


approaches  unity,  allowing  the  sum  to  diverge  as  more  and  more  terms  are  taken.  The 

2 

degeneracy  g  of  the  state  having  principle  quantum  number  n  1b  2n  [15].  A  method  for  avoiding 
this  difficulty  is  to  truncate  the  partition  function  sum.  The  amount  that  the  energy  level  of  the 
last  term  taken  in  the  truncated  sum  lies  below  the  ordinary  ionization  limit  corresponds  to 
the  lowering  of  the  ionization  potential. 

Several  criteria  have  been  put  forward  as  to  where  to  truncate  the  partition  function  and 
how  much  to  lower  the  ionization  potential  [4,  16,  17  ].  The  actual  value  of  the  partition  function 
is  rather  insensitive  to  the  particular  cutoff  criterion  used.  The  partition  function  sum  will  be 
truncated  when  the  classical  semimajor  axis  of  the  electron  orbit  is  greater  than  or  equal  to 
one-half  the  Debye  screening  distance.  The  results  obtained  here  by  physical  reasoning  agree 
with  the  Debye-Htickel  theory  for  dilute  plasmas  [18].  Good  experimental  confirmation  of  this 
Debye  distance -truncation  criterion  for  argon  plasmas  has  been  shown  [19]. 

The  semimajor  axis  of  a  classical  electron  orbit  is  [15] 


where  C  to  the  effective  nuclear  charge  (?  ■  1  for  atoms,  ?  *  2  for  singly  ionized  ions,  etc.), 
e  is  the  electron  charge,  and  W  is  the  energy  of  the  electron  state  in  ergs  measured  downward 
from  the  ionization  limit.  If  we  measure  energies  upward  from  the  ground  state,  which  lies 
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Wj  below  the  Ionization  limit,  we  replace  W  by  W.Qn  -  W.  Replacing  energies  W  measured  In 
ergs  by  term  values  E  measured  In  cm  *,  E  =  W/hc 


(14) 


Evaluating  constants  we  find 


a 


5.80699  x  10 


-4 


ion 


-  E 


cm 


(15) 


as  the  semimajor  axis  of  a  state  having  term  value  E  of  an  atom  with  ionization  term  value 
Elon' 

The  Debye  radius  is  given  by  [14,  18  j 


A 


d 


f- 


kT _ 

*£ -  D2 


6.90087 


/e+£"l(V1)2 


(18) 


since  we  are  considering  the  screening  of  a  stationary  charge.  The  sum  1b  over  all  the  positive- 

ion  species.  If  only  singly  charged  ions  are  present,  n  =  n 

©  + 


*d  =  4.87965  V'T JrT 


(17) 


If  we  denote  by  E^  the  term  value  of  a  state  having  classical  semimajor  axis  a  =  A^/2,  then 
from  Equation  15 


1.16140  x  10 


-3 


ion 


cm 


-1 


(18) 


Note  that  E  .  may  not  correspond  to  any  actual  term  value  of  the  particle.  Since  1cm  * 

-4 

1.23977  x  10  ev,  the  ionization  potential  lowering  is  [18] 

AV  =  1.23977  x  10'4(ElQn  -  Ed)  -  1,43>^  X  •  °  ?  C  ev  (19) 


For  plasmas  having  only  single -charged  ions  present, 

AV  -  2.95075  x  10‘8  VfT7TC  ev 
0 


(20) 
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where  ?  =  1  for  the  atoms  and  $  -  2  for  the  positive  Ions.  The  ionisation  potential  for  a  neutral 
atom  inside  the  plasma  is  therefore 


'ion  ’0 

The  Debye  -Hiickel  pressure  is  [18,  19} 


V,  *  Vn  -  2.95075  x  10'8  \TJ' f  ev 


kT 


24sA . 


(21) 


(22) 


For  single  ionization,  using  Equation  17, 


3/2 


-26  ne 

pd  =  -1.55517  x  10  "  atm 


(23) 


The  total  pressure  is 


p  -  nkT  -  kT-  « 
24tr*.a 

d 


(24) 


where  n  is  the  total  particle  density. 


PARTITION  FUNCTION  FOR  COPPIR  I  AND  II 

Knowledge  of  the  partition  function,  (Equation  3), 


z  ■£  >.* 
m 


■EJkT 


for  a  substance  enables  us  to  calculate  all  thermodynamic  properties  of  interest,  using  statisti¬ 
cal  thermodynamics  [5,  8,  20], 

Because  of  the  dependence  of  spectral-line  intensities  upon  the  populations  of  the  various 
energy  levels,  the  partition  function  enters  directly  (Equation  5)  into  the  calculation  of  these 
line  intensities. 

The  calculation  of  the  internal  partition  function  is  not  difficult  but  is  rather  tedious  for 
elevated  temperatures  where  large  numbers  of  energy  states  must  be  taken  into  account  and  an 
iterative  procedure  followed  for  an  accurate  calculation. 

2 

As  noted  in  Section  4,  the  degeneracy  of  a  state  with  principal  quantum  number  n  is  2n  . 
Examination  of  the  number  of  energy  levels  listed  by  Moore  [7]  shows  that  many  may  be  missing. 
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No  attempt  was  made  In  the  work  reported  here  to  use  lsoelectronlc  sequences  [7,  17,  21]  to 
add  them. 

Levels  above  the  ftrst  Ionization  limit  for  Cu  1,  at  62317.2  cm*1,  were  neglected.  Lines 
arising  from  states  above  this  limit  are  quite  broad  [22]  because  of  the  short  lifetimes. 
Spontaneous  lonizul  on  has  a  high  probability,  and  for  the  purposes  of  the  partition -function 
calculation  an  atom  In  these  states  was  considered  as  Ionized.  This  point  requires  further 
study. 

Figure  2  Is  a  plot  of  the  internal  partition  function  of  Cu  I  and  Cu  II;  the  data  are  taken 
from  the  appendix.  The  values  of  the  partition  functions  are  functions  of  temperature  alone  for 
the  range  of  conditions  covered  here. 


FIGURE  2.  INTERNAL  PARTITION  FUNCTION 
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SINGLE  IONIZATION  OF  AN  ELEMENT  IN  AN  INERT  DILUENT 

Consider  1  mole  of  element  M  undergoing  the  Ionization  reaction 

M  a  M+  +  e’ 

together  with  a  moles  of  an  inert  diluent  X.  The  reaction  is  written 

cvX  +  M  -  aX  +  aM  +  bM+  +  ce 

where  a,  b,  c  are  the  moles  of  M,  M+,  e*,  respectively,  at  equilibrium. 


(25) 


(26) 


We  assume  each  species  acts  as  a  perfect  gas,  with  partial  pressure  of  species  s  given 
by  [5] 


ps  =  V 

where  xg  is  the  mole  fraction  of  species  s  and  p  is  the  total  pressure . 
The  equilibrium  constant  is  written  [5] 

P  P  xx 

p  Pm  xm 


(27) 


(28) 


where  p  is  the  total  pressure  in  atmospheres,  due  to  particles. 
Let  N  *  original  mole  fraction  of  species  M 

N  m—1— 

1  +  a 

We  require  charge  neutrality,  so 

b  =  c 

The  total  number  of  moles  at  equllbrium  is 

a+a+b+c*M_ 


(29) 


(30) 


(31) 


Conserving  total  mass  of  species  M, 


M - aM  +  bM 


(32) 


Assuming  M  ■  M+,  l.e.,  the  atomic  mass  of  species  M  equals  the  tontc  mass  of  species  M+, 
Equation  32  becomes 

1  ■  a  +  b  (33) 
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The  total  number  of  moles,  Equation  31,  becomes,  from  Equations  29,  30,  31,  and  33, 


MT=N+b 


ai  equilibrium.  The  neutral  mole  fraction  is 


and  the  ionized  mole  fraction  is 


a  1  -  b 

Td "  mt = rm + b 


b  b 

X+  =  mt  "  (l/N)  7  b 


(34) 


(35) 


(36) 


x  =  x 
+  e 


by  Equation  30.  Let  x  *  x  ■  x,  and  so'ving  Equation  36  for  b, 

4*  6 


b  = 


N(l  -xi 


so  that  Equation  35  becomes 


(37) 


(38) 


XM  ■  1  tl|1,~ ’‘i  ■  M  -  (N  . 


1)X 


N  N(1  -  x) 

Substituting  Equation  39  Into  Equation  28  yields 

2 

v  x  p 

p  N  -  (N  +  1  )x 

Solving  for  x  and  taking  the  positive  square  root  since  x  >  0,  we  find 


K  (N  ♦  1) 
x=-£ 


Tp- 


’/♦- ^-,-1 
y  k  (n  ♦  ir 


Whenever 


Equation  40  can  be  written 


c - SEE— «i 

k  (n  + 1  r 

p 


J£Lri.Ic+-i-c2+  1 

N  +  l|_2  8  C  +  16  C  +  '  ’  J 


(39) 


(40) 


(41) 


(42) 


12 


Institute  of  Science  and  Technology 


The  University  of  Michigan 


The  neutral  mole  fraction  Is  then  given  by  Equation  39. 

The  inert  mole  fraction  Is  obtained  using 

x  =  1  -  (x  +  x  +  xM)  (43) 

a  +  e  M 

To  obtain  we  proceed  as  follows.  The  Saha  equation  (Equation  7)  becomes 

n  n  (2jrm  kT)3//2  2Z  -V./kT 
+  e  e  +.0 

— —  - * - 5 — e  (44) 

nM  h3  Zm 

where  VQ  is  the  lowered  ionization  potential  of  the  neutral  atom,  n+,  ne,  nM  are  the  number 
densities  of  ions,  electrons,  and  neutral  atoms,  and  Z+,  ZM  are  the  Internal  partition  functions 
of  ions  and  neutral  atoms. 

Using  the  perfect  gas  law  for  the  partial  pressure  of  species  s  [5], 

ps  -  ngkT  (45) 

We  can  rewrite  Equation  44  as 


P+pe  3 

—  «  (2irm  ) 
PM  ® 


-vkT 


3  z. 


dynes  cm 


to  have  the  same  form  as  Equation  28.  Changing  pressure  units  in  Equation  46  from  dynes  cm 
to  atmospheres,  and  evaluating  constants  gives 

.  /22Z  -11605.4V /T 

K  -  (3.28978  x  10  7)T8/Z~  e  0  (47) 

P  ZM 

where  Kp  has  units  of  atmosphere,  T  Is  measured  in  °K  and  Vq  Is  measured  In  electron  volts. 
We  can  also  write  Equation  47  in  the  form 

/2Z  \  VQ 

log  K  -  -8.48284  +  }  log  T  +  log  —  -  5040.16-=r  (48) 


The  total  pressure  Is  given  by  Equation  24: 


p  ■  nkT  - 


Noting  that  n  ■  x  n,  and  evaluating  constants,  we  can  rewrite  Equation  24: 
€  © 


0.734002  x  1022  •!+  1.14150  x  10 


-4/  e 


-  n  -  0 
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Since  the  Debye -HUckel  correction  must  be  small  to  be  valid,  an  approximate  root  of  Equa¬ 
tion  49  is  given  by 


n  =  p/kT  =  (0.734002  x  1022)p/T 


(50) 


This  estimate  may  be  improved  using  Newton's  method.  Thus 

[(0.734002  x  1022)p/T]  4-  [l.!415  x  lO'^x^/T)3^2]  - 


V =  nj 


(1.71225  x  lO'^/T)3'^  -  1 


(51) 


where  n^+1  ie  the  improved  value  for  the  root  n  of  Equation  46,  using  the  results  of  the  j-th 
approximation,  n^. 

After  satisfactory  convergence  of  Equation  51  has  given  a  value  of  n,  the  electron  density  is 
given  by 

n  =  x  n  (52) 

e  e 


from  which  the  Debye  radius  and  tonizatton-potential  lowering  Equation  19  can  be  computed, 
since  the  temperature  is  known. 


7 

APPLICATION  OF  SPECTRAL-LINE  INTENSITY  MEASUREMENT 
TO  TEMPERATURE  AND  COMPOSITION  DETERMINATION 

As  a  simple  illustration  we  consider  the  single  ionisation  of  an  easily  ionised  impurity  in 
an  inert  carrier  gas.  We  shall  measure  an  impurity  atom  to  impurity  atom  line  ratio  Rq  and  an 
Impurity  atom  to  impurity  Ion  line  ratio  R+.  The  subscript  zero  denotes  the  neutral  state,  and 
the  subscript  plus  the  ionized  state.  Let  the  two  neutral  atomic  lines  arise  from  transitions 
with  upper  upper  levels  Eml  and  Em2;  the  tonic  line  has  upper  level  Em+.  We  shall  assume  the 
values  of  u,  gmA™,  and  Em  are  known  and  that  Z^(T),  Z+(T)  are  known  functions  of  the  tempera¬ 
ture  T. 


From  Equation  5:  the  line-intensity  ratios  are 


and 


R 

+ 


exp 


kT 


-(E 


m+ 


kT 


+ 


(53) 


(54) 
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where  Vq  Is  the  Impurity  ionization  energy.  Note  that  Em+  Is  measured  with  respect  to  the 
ground  state  of  the  Ion  which  lies  Vq  above  the  ground  state  of  the  atom.  If 

(hm  k)3^  . . 

a  - - - =  2.41469  x  101  (55) 

h3 

then  from  the  Saha  equation  (Equation  7)  and  Equation  8, 

JL.f!l0T3/VVkT  1“' 

no  h 

Equations  53,  54,  and  56  constitute  three  equations  in  the  three  unknowns,  n^,  n+  and  T.  If 


and  Equation  53  is  solved  for  T, 


, .  ‘'■(vC)i 


H  -  z/z0 


E  ,  -  E  . 
ml  m2 


k  fn  (a0/R0) 


_  n0Sfl  -<EmrEm+-V/kT 

R+ 


n;  3/2  -V0/hT 

—  -  2a0T3/  e  0 
n0 


From  Equations  61  and  62,  we  find  the  Impurity  ion  density  to  be 


Equation  54  then  becomes 


and  Equation  56  becomes 


2aR+T  -<Em*-Eml+2V0,/kT 
n  - - e 

+  »♦ 
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In  viewof  the  charge  neutrality 

The  neutral  Impurity  density  Is 


n 

+ 


=  n 


e 


2  VkT 

n  e 

+  _ 

2a(3T3//2 


(64) 


(85) 


A  further  measurement  of  a  carrier  gas-impurity  atom  line  ratio  R  would  give 

c 


„  'c  ncz0 


(66) 


where  the  subscript  c  =  carrier  and 


"I  (*mA”)l 


(87) 


Since  T  Is  known  from  Equation  60, 


n  *  ■ 1  'L  e 
c  *cZ0 


The  pressure  is 


(68) 


P  "  ("c  +  n0  *  n+  +  ne^kT 


(69) 


The  mole  fraction  of  neutral  Impurity  is 


0  n„  +  n  +n  +  n 
0  +  c  e 


(70) 


and  the  mole  fraction  of  Impurity  ions  Is 


♦  n„  +  n+n+n  e 
0  +  c  e 


(71) 


The  total  mole  fraction  of  impurity  is  x+  +  xQ. 

Note  that  the  assumption  of  thermal  equilibrium  permitted  the  same  temperature  to  be 

used  In  all  equations.  Note  also  that  E  and  E  ,  must  be  referred  to  a  common  zero. 

me  ml 
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S 

METHOD  OF  CALCULATION 

Equations  3,  16,  20,  40,  48,  and  51  were  programmed  for  an  IBM  7090  In  MAD  (23),  accord¬ 
ing  to  the  following  scheme.  The  pressure,  temperature,  and  original  mole  fraction  were  read. 

The  partition  functions  for  Cu  I  and  Cu  II  were  calculated  by  truncating  the  partition-function 

■6 

sum  when  the  terms  became  less  than  10  .  From  these  trial  partition  functions,  a  trial  K 

P 

was  calculated,  assuming  no  ionization  potential  lowering.  This  trial  permitted  a  preliminary 
calculation  of  ne,  Ed,  and  These  values  were  used  to  re-evaluate  thepartltion  functions,  and 
the  sum  was  truncated  when  the  terms  became  less  than  10'6  or  the  term  values  were  higher 
than  the  lowered  ionization  potential.  Second  values  of  Kp,  ne>  Ad,  and  lowered  ionization  poten¬ 
tial  were  then  obtained.  When  successive  Iteration  of  this  process  resulted  in  a  fractional  change 
-5 

in  Ad  of  less  than  10  ,  satisfactory  convergence  was  assumed  and  the  rest  of  the  particle  den¬ 
sities  and  spectral-line  Intensities  were  calculated.  New  values  of  pressure,  temperature,  and 
mole  fractions  were  read,  and  the  process  was  repeated. 
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Appendix 

COMPUTER  PROGRAM  AND  TABULATION 

Values  of  the  Internal  partition  function  (Equation  3)  for  neutral  copper  (Cu  I)  and  singly 
ionized  copper  (Cu  II),  the  ionization-potential  lowering  of  Cu  I,  Debye  radius  (Equation  17), 
Debye  pressure  correlation  (Equation  22),  equilibrium  constant  (Equation  47),  particle-number 
densities,  and  absolute  spectral-line  intensities  for  Cu  I  >.  5106  and  Cu  I  K  5153  are  presented 
as  functions  of  temperature,  original  mole  fraction,  and  total  pressure. 

2000  <  T  <  10,000°K 

0.03  <  p  <  10  atm 

0.0001  c  N  <  0.03 

The  number  .2002  x  101  is  printed  out  as  .2002E01,  and  .1192  x  10*1  is  printed  out  as 
.1192E-01. 
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Alls,  ilii'irtii.  fMtJS 

CP. AS  NO.  IS/ 

01  NANAI  NCI  UBS  COAACNAI I  UN 

oi msm  *y»iia»  nivliirv 
•ui  1 

IASI  A  MANMANA,  LAI  IM1AS|A 
AMS.  ACAl.  NUNINl  I. 

I  IH4A4 | AN 

CcHY  NO.  In  1 


mis!  IS  A | NLNSr I  Cl  HAAN* 

shah  mirn  oivimun 

A  a  C.  Hfl  veils 

Ml'.  ASUIllSi  LAL  II I1NNI  A  SlOl  V 

AMS.  SlHOUM,  *10-4*0  t  . 

CI.AV  NU.  HI 

AahCUH  MIVAMCh  IlCSIlAMOS 
AhVkllt  til  VISION 
U>  B»SI  I’.lh  HNtl  I 
CM  I C  At  II  is,  ill  INCH 
AIIS.  ('OCUMINI  l  |l  MANY 

CUAV  Nil.  /CV 

i'.si  l lull  rra  nirisLi  analysis 
IbbA  rilSSICIIUM  AVtNUli  >*.  A. 

NASI-  I  SC  I  IIS  V.  0,  '  . 

aiis.  uiiiinnaNi  mu  is  n, 

f  DAY  NO,  VIC 

IN. 1  Hull  YOb  Pt  AIN 11  ANALYSIS 
Itibh  LOSNICIICUI  AvINUt.  N.  b. 
bA  SI- 1  SC  ION.  r.  c.  /coot 
AIIS.  CNiSllM,  kAHL  M. 


UllCUNINI  LOS  I NOL 

IN, I  Hull  ION  OlAISM  ASALVStS 

ISbB  CI'AAilI iLii  t/flul  N.  s. 

u An  t suns  v,  r .  i, 

ills.  IttiLhk  H  A  i  ON,  A  ('CL  A  N. 

LOW  Nil.  / 1/ 

ISStllUll  MIN  PA  V I  SSI  ANAL  YUS 
Lb*«  CIlSNI  Cl  |  (.0 1  AVI.,  S.  N. 
NASI'  I  SI.  I  UN  1,  II.  I  . 

AIIS.  1 1  Liu  |SA0N/AI|(lS  01  A  JU 
Cl  NY  SO.  /II 


Cl  St  A  *1  UICIMIC  IflPAASV 

I  CMC  I  IVIMUS 

ns  hah  slant 

S  AS  I  A  HANNAN  A,  CAL  |A||ASIA 
AIIS.  UUCuMISI  LOS  I A  01  OtflllM 
CLAY  NO,  ns 

Qt  si  aii  me  Ink  i  oaaasy 

I I  A AC  UlVIL ION 

ns  1 1 A  1 1  St.,  A.  I),  IINAblN  VS 
S  AH  I  A  NAHBANAt  CAI  IAIINMA 
Alls.  I  ICAMCm,  At  Ak  S 

UlBlCfO**,  UASA  UAH  (IN. 

CCAV  NO.  I H 

r.lSIAtl  met  MIC  HiPAANV 
tl it  tans ICt  rmsK  s 

SY>ACUt  I  Nik  Y  C  MM 

Alls,  i  Inn An  I  AN 

LLTY  Nl),  IH 

lit  VIA  M  mtl-U  I  OAAASY 
L  Ibht  A  ll  II  AAY  ILH.IS.  OtAI, 

AL  I  (iOAIf  Mill 
ill  1C  A,  Stb  »«H« 

Alls.  LlbiAAY 

loav  Nil,  Hr 

A.  L.  AC»  IAA 

A  |  SL  CA  AAUttlA,  MSNLYi  VANIA 
AIlN,  A  INL,  0«.  b.  L, 

LCA»  NO.  Its 

Al >101*1  ANVtICt  I  IANSNT 
Cl  Si  AAk  ILICtllC  LOAAASY 
•  lltlll  «  IAALI  VthILkt  MAI, 

*AAtl  ttlhNCHISV  MAN  ANY  i».  |IA»I 
VAlllY  AIM  lit ,  AlNNSHVANIA 
AMS,  |  |  A*  AH  I  AS 

LOAV  NO,  l»b 

wlSlAAl  UHIAI;  LOAAASY 
Al  rstAV  tvtltbi  tret. 

IIIN  CHfttNUl  11*111 
Ah  ll AtriAhlA,  AA.  I  MCA 
III*.  ACI4V  ,  JAtOAul  A. 

re  ay  no.  itl 


Cl  Nl  A  A|  Ufcllll  II.APASY 

«A-/str*  mu«i  iiai. 

Ill*  Ihl-'-Ngt  '.lain 
Ah  |l  ACIL  Aha  *  AtNSkUVANIA 
AIO,  •  *  CONICAL  ISAONAAMIlN  LASIIN 
CCAY  Nil.  I  A.* 

CISINM  ILICtN  1C  CCAAASV 
•ISSIU  A  1AACI  01 Y I LICS 
SAACI  ftCh,  CIA  A  a  U.  Nfll  HS*S 
ANIL  AltLAh  14  I,  AISNAYi  VANIA 
Alls.  rCCUOtNt  llAMANY,  AOOA  l|l«l 
CA  All s.  i Ann  t  SC  f  I. 

COAV  NO,  la  I 

OINtMl  menu  CWAANY 
■  mm  a  AAAct  o i v i lies 
SAACI  fICh.  CfNIlA  A .0,  AO*  »»S* 
ANIL  AtllAl  IA,  AA.  LVIOI 
AIIN.  AAOblN,  JOWAA 
ICO*  ■* AM 

COAk  NO,  INA 


CtSIAAL  met  A  1C  COAAASV 
■lYL  III  A  tAACt  OIVUICN 
LA4CI  IKK.  CANltA  A.Q.  NO*  NISI 
ANIL ACIIAAIA  ,|.  AANNIHVANIA 
AHN.  Mill.  AAteiCN  4. 

COAV  NO.  le* 


Cl  SINAI  AL  IONS  ClWHUbA  1 1  ON 
OHISLI  NISfAMLA  I  AHOAAIIIH  1 1  S 
»U«  1 

SASIA  ■  Am  MAM  A  ,  LAI  |tllNS|A 
AIIS.  UCUl  i  mHm.  41  i. 

Cl, AY  NO,  IS* 

lit  SlN  AL  Al  IPNS  CIINAL'NAIII'N 
III  A  (  SSI  NttlANI.h  I  AHUM  A  tl>B  II  S 
NIK  I 

SASIA  HAmmAmA,  LAI  IIIIHSIA 
Alls,  AAltrs.  LLll'i  J. 

I  I  AY  Nil.  IN'. 

Cl  St  SAL  ANIUMCN.  I  SCI  NAlWAtl  * 

1 1  anASCOAI  flvislis 
At  m  at* 'INN  AVINUI 
ClINfAII  t.  CAl  lAi'NSl A 
AMS,  SAMS,  SAtllAS  4, 

SLAINV I  VIA  l  |  HA  Am  t  Slav. 

cr.AV  NO,  |1» 

l  |  hN  AN  |  AS 

til  SI  NA|  AMILIStCN,  IMCaAONAII" 

M  ANA  C I  I  rivISICS 
l/.S  ALVan  AVlNk" 
l  H  ll  I  »AllS,  stb  JlNMY 
SMS.  HAM.btN,  ASAV  S, 

U’AY  SO.  IS  I 

lit  uAhYS  US  LONACbA  I  (UN  III  AA.mIIA 
NIIAL  IStiltlS  NOA| 
bicrcir,  mass AY/u. im 
Alls.  I  •  IhNIt  At  L  |  Hi  ANY 
SA/IAAD,  L. 

CCAY  NU.  ISA 

0( UAhVS ICS  COAAUNAIION  Ot  AA| N |(A 
ANY  IKS  Nl  ll  AMth  I  lYlllllS 
•UNL  IStiltlS  AIIAI. 

AlCACHO.  bAUACAUll  IIS 
Alls,  AAktlAAS.  JINCNt 


bi assini  instant*  lcaacaaius 

ileeee* 

Ib.o  tn,  AOUNI A  IN  AVtNir 
nuAAit.  caliaonnia 
AIIN,  LuNACt,  CICIlY 

CCAY  NO,  4lO 

CUUCVfAN  AIACAAM  CCAACA A 1 1  IS 
L/I*  Alls  (ICS  N04I 
aancn  I*.  OA  in 
AIIN,  I  SO  INK  N  ISO  llesAAV 
P/INA-O 

COAY  NO.  it  I 

umuanas  Aiacasn  t ShInacn |Nt.  lonh, 
ay  teste*  rtvmsN 
bl I AAACI,  Nib  vlnv 
AIIN.  •  SCtlNIi*  ISO  L  INN  ANY 
AlANt  * 

CCAY  NO.  iC» 


MAY  Y  S  AINCNAM  CON  AC  A  A  I  ION 
AO* t  CM  ICt  Ml  Mtt 
■  |4A|S0»AS  A,  At  A*  AAA 
AIIN.  CNOa.  NON  I If  A 

(AC..  IICMNICAt  UN. 

COAV  NO.  lb* 

MILICCVAI  CCNAQAAI ICN 

its*  bittbOCC  eCMCVAer 

10*  ASOtill,  CALIAONNIA  NOOA* 

At  IN.  A  lies  AS,  c«,  IAUL 
ttcueitv  oar  kin 

CCAV  NO.  IOA 


Alb!  A  At  Cl  AC  I  H|<  CUb  A  AAV 
■mill  *  * AAC I  01  VlllCN 
tAAC!  Hth.  Ctstlb  A.  C.  10*  1*1* 
AHlICflAhlA,  AA.  Midi 
AttN.  OUtf PAN,  AM  A( a 

COAV  NO.  INA 

AtNfAAi  ILACIAIC  C0AA4SV 

smut  •  *aac i  oi vt i its 

tAACC  MCA.  ((*»•  A.e.  to*  Alt* 
Ahll AtILAh IA  I,  ArSNtVL VAN) A 
Alts.  It*.  Afsav  I.. 

•A.  «IAK 

COAV  NO.  t*t 

WMIH  mete  1C  COAAANV 

■miur  •  »AAct  oi  vines 

IAnCI  MCA.  CINMA  A.  C.  BO*  «**« 
AhllACIlAUA,  AA,  IM«| 

AttN.  A IIIMQ* ,  I  NBA A*  t. 

VIA.  al-CStNV  A  MV  1 1CL  IIAIAIV 

coev  so,  tee 


Mee-tltSIS  ISiONAli  A  A  1 1  c 
A.  c.  SCI  A9 

*IAM  CClltCC.  AISSIVLVANIA 
AttN.  IIANANV 

COAV  NO.  tr \ 

UtANICAL  flC'-A/SM  List l« 

♦NIL*  I  »  A|NCNA."  Cl-hAANV 

*-a;c  mm*  ot  via  ion 

tUilMI  AVlSuf  »  HAL  I  ttaett 
CUlVie  CltV,  CAL  IVCNNt  A 
AttN.  Ullllllt!  •.  C. 
HAfavitO* 

nth.  UOCUMISI*  emu 

CCAV  NO.  ICA 

MU CAI ft  AINCSAtf  COAAANV 
MOUNC  tVIMN*  000 UA 
AWU.ietCN,  CALIIOANIA 
•  TV'.-  MA|B( .  A  AAV  ICU 
HCMNlCAt  1 1 Naaav 

COAV  NO.  >01 


OULU At S I  COSIHUL 
<  KSli  lull  AON  lltMSVl  ANAL  VMS 
-  LAN  CliSNfUIU/l  svisur.  S.  k. 
bA>l  ISutllS,  P.  C.  ii'vOH 
Alls.  aahI is,  JUHS  4. 

VIA.  Dl.CuMIN  I  COS  I  Mill 

LOAV  SO,  >1* 


niiLUAtsi  cos i niil 

I', Sl  I  lull  Alia  MAI  SSI  ANAL  *  Si  S 
|  Hblt  LLSllCtllUl  AVINUI.  N.  h. 
NHSh  I  S(.  I  US,  I  .  L.  >»C'/S 
aiis,  ■  •  s» r s ,  Jutt-LA 
vis,  nc’Luaisi  i.Oshui 

LOAV  NO,  ilk 

|  .Nl  IIUII  Ilia  PAH  SSI  »NAl  V SI  S 
K.LW  CrSSICIIlUl  AVINUI,  N.  a. 
NASA  ISl.lUS,  I,.  L.  iCLLS 

ATTN.  eCKULTU  WILLIAM  J. 

via.  muAisi  cos i hih 

•  I.AY  NU.  /IS 


I  .s H tint  nr  "■  * v ai  slut  it* 

S*.  If  I  hSnlllLY  tei'ANI 
cat  AA  l<  (I  I*.  "HtSAll'Uk!  Ilk 
*• ih .  I  IHHANY 

LOAV  NII,  /H 

1 1 1  H  ClIH  AIINA I  ms 
VLvVA  rtVltlUH 
I*. 4  Aid  AlLl  « HA  I 
AAlll  At  IP,  LAI  I A  04  S  I A 
At  I  4  ■  OLUlMtSl  ClA,  1 4111 

Lt'AV  r  /la 

I  III  L114AHNAI  ICS 
V  li-f  A  ll'vltus 
l*w  aali  am  HUAI 
A  alt'  A|  111,  CAI  I  f  Oh  S  I A 
Alls.  milUAisi  costaui 
4ALINI,  AO*  IS 

OI.AV  so.  ill 


lit  a  i  pi  a  Ai  i  AaCBAioair* 

AL.  NASh  ISti II.S  AVI  NL( 

Null  I  V  lb,  SIB  4f*  SI  V 
Alls.  L  IahANY 

C PAV  NO.  m 


ilil'S*  ht'ANIS*  usivtani* 

AAAI  lie  mil's  LAACNAIONV 

a*/ 1  CtiadA  avcnui 
klivca  laaisc.  AA*viAsr  |w«u- 
Alia.  tea  IN,  OA.  /Ann  A..  JN. 
VIA,  nCCUMINIt  UANAAV 

COAV  NO.  »l 


40AM*  MCAalb*  UNlvfantv 
AAAI  110  PAVnt*  LAAOAAIOav 
esil  CIcaUA  AVINUI 
*  lives  SA4ISCV,  MA4YLAS0 
AIIS,  L  lANAa |AN 

COAV  NO. 


in 


■  MAAS  AlaCAAtt  C0NACNA1  ion 
NULL  * *•  CIVH10N 
OAACIS  Ct  Iht  600*  SOAL 
COiCAAK  *AA|NUI,  LCtOaAOO 

aiin.  ciAtllAlft  rutgntst  eaiscn 
LOAV  NO.  lit 


l l AN,  ISCOBAOAAtei 
IH  ICNI A  AVINUI,  N.  b. 
ONASC  NAAtOt  I,  AlChlOAN 
AVIS.  *|At,  LUIS  A. 
u SHAN  I AN 


CCAV  NO. 


II* 


lilies  mllMt,  INCCNACAAtfC 
III  cat  lust  klatt  I 
MAltAAM,  NALtACAUtMl* 

Alts.  *te(*IVfN,  CtCNCf  s. 

CCAV  NO.  1,51 


LOCAMf  IC*(Al  |  IOBN I  A  CUAAASV 
fi-i*/»lLL.  Al/AIANI  A*  I 
ewaeANa,  calKonnia 
At  IS.  CISISAl  limasv 

aUfftCMM,  ON,  tCBAA  b. 

COAV  NO.  lib 


inClHfll  CA*  IVOBNI A  COAAANV 
A.  0.  *CI  e»l 
twACASa,  CAiirOBMIA 
Alts.  null.  AAA  ION  l. 

CVKf  l  ISAM  IAN 

COAV  NO.  li 


lOlaHUl-  MACNAM  LCNAr*nil,N 
■  mill  •  SHALL  UIViltCN 
hAsnvta  tint  i 

All  C  Ulti,  CAllfOANI* 

Alls,  aC/UAAlU,  I  a,  NA.  A,  Him 
KGH.  IShCNh-IIOS  LfNItS 
CklVV  NO.  .It 

t  ilLKHMl  HAS  Ilf  •  SVACf  COAAANV 
A,  C.  ikCi  SO* 

IWSSVVAK.  CAL  IfCSSIA  SA'JAA 
AIIS.  PAH,  AUtCS 
nl  *I-Al 
VIA.  Nlft-NU 

CklAV  sn.  141 


U  AN  ASIAN 

LUCaMtll  mCHOSILS  LfMAASV 
U,  I.  MKhbAV  it 
AlAINUIir.  Sib  JfNtlV 
AHS.  CbtlNUSC,  H  lltkAAp  C. 

COW  NO.  ill 


OllLUAINI  cnSUIIL 

UJIAL  ILICMONKk  CCaAPNAf ICN 

a/*  vaiisk  alvtM  isisut 
Nib  vc aa  ti,  vi«  vnac 
Alls.  SIAL'N,  ALLIS  I  . 

CHAV  Nfu  4 II 

ia«»H  llllh,  altlAHCM  INtlllLH 
LllUILL,  HilAiKtU  III 
ills.  NieasnA.  a.  a, 

Cl  AV  NO.  *  li 

■axils  CUAAASV 

oisvia  tivinns 

A,  Ci  Nf.k  Ills 
ftisvfa  I,  CIILIINADL) 

Alls.  ALVASCID  KChSUirtiV  L IhAAbv 
LOAV  111,  ill 


MAN  I  IS  CUAAASV 

aos i  pm  ici  eni  hi 

OtSVfN  I,  CCLL4AIKI 
AHS.  Ml  COMNlia,  4*t»  A.  ACS  LSI  I 
PA  I  AC  N(S.  t  IA.  !•*/ 

CCAV  NU.  ilk 


t,  IhHAm  |  AS 
MAM  I  IS  CI.MAANV 
Ansi  CM  ICI  AUl  tbit 
UNlASt  U,  ALCi  I UA 
AIIS.  liAhLUS,  APSA  C. 

CHAV  NU,  i  AS 

■«Hl|S  HAM 1 1 M  A  LU4AU*At|UN 

Iai 1 1 ■nai  i,  MANYiasu 

Alls.  SUINCI  fLChkUOCV  1  I  A- AMY 
CCAV  NO.  Ml 

NAkkACUI.*lltt  IS*I.  Ilf  IlLMSdlCCY 
|  ISCUIN  L  ANObaitlNV 
i*«  bC.U  MAftl 
|  k  »  (SI, IDS.  AAUACUUlf  HSnilM 
AHS.  HAk  SM,  MAat  IS 

COAV  Nl),  /It 

■astAcuciMt*  in*i  .  nr  trCNSuicuv 
I  I  seas  l  A»OMAICHY 
kinSCICN  f)|  H«*taCMU*im 
ails.  ■luuAt.  *lb> f N I  4, 

CUAV  NO,  MS 


■anacnckM is  ism.  «»  mcnsiilixv 
i  i seas  laer-HAiObY 
aiis  I  mini  Ana  n 

L,  H  INL  ICS  M.  •maCHUMIH 
ai  is.  Ct  ana  ,  rnsAi  r  i. 

CCAV  NU.  Ilk 

■assacHckt 1 1*  ism.  <>h  i|lnsiiii;iiv 
L ISCOLN  I AbONaluev 
A.  C.  be*  ft 

1 1  ■  INCtCN  tl,  AAllAlMkiStm 
aiis.  fuueiac,  \i  vaou* 

CCAV  NO.  i*b 


HAktachun  1 1  *  Ism.  uf  lAiMsui  rev 
I  INLOLN  LAblliAlUav 
au*i  6»f it#  eo*  ti 

LfllSCfCS,  AAkkACICklft* 

AHS.  AUNItli  AANV  A, 

IICCUAINM  L  |  »B  A*  I  AN 

CUAV  NU.  i* I 

AA**a(MC*if tl  ism.  uf  itCHsu.'-c* 
LINCOLN  LAbCAAioav 
A.  0.  aca  ft 

IIAINCION  ft.  AltkACMUtf  ML 
AttN.  Af  SAAAAA,  fA AW  L. 

COAV  NO.  i«> 


NAIftACMUf  Ml  ISM.  Of  tfCMSIliliav 
UNCO* 4  l Asoeitoev 
Auftt  C»Mt«  eo*  M 
unsetcN  ti,  NAMACMuicm 
•ITS.  AlAAtet,  C*.  6*1*  f. 

coev  no,  is* 

N*»»*Cl*»tMft  INM.  UA  tfCMSOi DO V 
LINCOLN  IIK'UQM 

is*  veot  ftfeetT 

if*  IN6T0N,  AAklAClClIItt  4IHI 
ATTN.  UATTtS*.  elCAABf  C. 

CCAV  NO,  i*A 


NCPONNf it  AiecBAfi  coe*o*Aiirs 

so* i  off ici  eo*  *i a 

It.  ICt/lt,  AlkftOWAl 
AHN.  IS4INtf*tM  LIBNAAV  6(Af. 
NAUM,  J,  a,  U|t! 

COAV  NO.  #** 

NtiAte,  iNCCAAoeetce 
toee  as  UNO  t  on  eouuvA*o 

A  ALL*  CMieCA.  VIAOISIA  IM*A 
AttN.  MCMItCAl  ISAO.  CINfl* 

coev  no.  m 


IlkMimit'MMNilU  ait.  CO. 
OUffMV  MOOUCtl  CONACtAtlON 
MMt*  CIV  MI  * 0*61101  fine 
Ml  ICN  It,  M1IACMi«6«n 
Atts.  t level tt a*  Lieseav 

COAV  NO,  lAt 
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•  1 1 ,  incuv*c**wq 
>iu  avinui 

c *m«h  1  oci  »*.  a*i,*tKim  m 
lllk.  UHtH 


m k  i ».•  *v .  inco»ao**ic? 

■icacaivt  ♦  1m**«*i  baciacign 
itu  auiaan  avpnw 
c*a#aipet  tv.  M|tACnut*m 
*1  IN.  IC**M.  I  AMlNCt 

CCAV  Nil.  1*1 


KUMiUt  MIIMCH  CUBAttAAIICN 

mi  aitiaooD  bouwcvmc 
wot  ihuii  I*.  c*upo*nia 
alia,  paavhnip,  uiiabiih 
i  laaaaUN 

CU'V  NII.  14* 

Pitas  law  aniaatn  ccaagbamcn 
nil  aftlartcr  aOi*  uttC 

wot  IMUh  I*.  CAUfO*NIA 
at  IN.  PC4N*  i  .  Nll)»ANN 

occuMCNt  coniaow 

CCf*  NO.  I»» 


WINN**!  Aiavlitt 
M I  r at  ica*0«AU0N 
NON  I  CNN  Hi  ao*  I(  * 
aiotctn.  N»tt*wfutrm 
Alta,  uarw  in.  jf*s  * . 
w tea** (an 

HIM  HU.  INC 


attat  tc*N(ta*i  ion 
•tCfCan,  AAttACxutl tit 
NUN.  I  lea*** 

MCNAUQA' >  10*4  t  .  IWtal 
UIAV  NO.  INI 


nq*ia  am  at  can  *vl*M04 
*01  cat  t  ICt  CIVltll'4 
N.  c.  ao*  a- i 
anaaUh.  caupovni* 

i»ir*ati.  »■»*. 


)****« 


:-•»*  no.  i»; 


no* it*  *ncmiw*n  ayianuN,  Ini. 
aowatmvNt  CIvltli'N 
N.  C.  ac*  «N| 

C4MCW*  **•«.»  C*gltS.SNI*  -II  tc* 

*1  IN.  DlatlW.  AAlVU*  0. 

CO**  NO.  INI 


i.  a.  iimn  AttuciAitt,  inc ■ 
1**4  NiitHiat  bOuuvaaw 
tANlt  *taic*.  gaimobnia 
at  in.  ttcuait*  ui*«i»n 

COA'f  NO.  IkA 

c.  ►  >  am  mu  *noti»ut.  inc. 
M«  aiwt*>l«  aiiuityaac 
t*Ni*  ncnic*.  catitoaNi* 
ana.  attNaai .  sabqwo 

C0*«  NO.  I*N 

■  4i. io  cc»au*amon  o»  A*r»it* 

*i net* *ct  con .  «  ccNiar.it  civ. 
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